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1 Introduction 



The recent evidence for CP violation (CPV) in singly Cabibbo suppressed D° decays 
to two-body final states from LHCb [1] and CDF [2] has heightened interest from 
theoreticians in charm physics. In order to investigate these promising hints further, 
multi-body decays with three or four particles in the final state will be needed. Direct 
CPV arises when decays with two possible routes to the same final state with different 
relative weak and strong phases interfere. The intermediate resonances in multi- 
body final states have different strong phases, and thus if there are also different 
weak phases, CPV is guaranteed. Moreover, in multi-body D decays it is possible in 
principle to extract information about these strong and weak phases which remains 
hidden in two-body systems. 

Furthermore, some of the easiest charm decays to reconstruct experimentally have 
three or four charged particles in the final state. For example, the most abundant 
flavour-tagged singly Cabibbo suppressed final state available to most experiments is 
the three-body D + — > K~K + n + decay [3], and so, naively, this should be the most 
promising charm decay mode for CPV searches. 

However, multi-body charm decays present several key challenges. To exploit the 
phase information to the full requires an amplitude model. Such models carry with 
them systematic uncertainties that are difficult to quantify and depend on choices 
made by the analysts. The properties of the many possible scalar resonances re- 
quired by the amplitude models are often poorly known, and model uncertainties 
now approach or even exceed statistical uncertainties in many amplitude analyses 
(e.g. JJ|). Therefore it is worth investigating model-independent techniques which 
are not subject to these uncertainties. With the large data samples now under study 
at the B-factories and LHCb, such techniques should also be faster and simpler to 
implement. 

There are now several model-independent techniques in the experimentalists' ar- 
moury. I will discuss recent analyses which illustrate some of these well. 

2 D° hhn° at BaBar 

One of the first model-independent searches for CPV was performed at BaBar on the 
.D* + -tagged decays D° — > 7r + 7r~7r° and D° — > K + K~7r° [5]. An amplitude analysis, a 
phase-space-integrated search for CPV, and two model-independent techniques were 
employed. In 82,468 ± 321 tt+tt-tt and 11,278 ± 110 K + K~ti decays, no evidence 
of CPV is found. 

In the first search, the Dalitz plots of these decays are divided into simple grids 
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of bins. In each bin, the variable Sep, defined as 



is calculated. Here, N l D+ represents the number of D + decays observed in a given 
bin i, AN l D+ its uncertainty, and a is an overall normalisation used to cancel any 
production asymmetry effect. This is similar to the "Miranda" method [5J. 

The sum of the squares of the S l CP over all bins is a x 2 f° r consistency with no 
CPV, with a number of degrees of freedom equal to the number of bins minus one 
(due to the overall normalisation). BaBar obtain a one-sided Gaussian confidence 
level for consistency with no CPV of 32.8% for 7r + 7r~7r° and 16.6% for K + K~n . The 
distribution of S l CP values over the Dalitz plots is shown in Figure [TJ 




Figure 1: Distribution of S l CP values over the Dalitz plots: T>° — > 7r + 7r 7r° on the left, 
V° ->■ K + K-it° on the right. 

BaBar also perform a complementary search for CPV which makes use of the 
angular moments of the helicity angle Oh- This is defined for decays of the type 
D — > r(AB)C as the angle between the momentum of A in the AB rest frame and 
the direction opposite to the D momentum in that same frame. The procedure is to 
bin the Dalitz plot projections in two-body invariant mass, then weight the efficiency- 
corrected number of events in a bin by the I th Legendre polynomial function of cos 6h- 
This yields moments 

p o= \S\ 2 + \P\ 2 ^ p i = v /2|s||p| cos ^ p) p 2 = y|| P | 2) _ (2) 
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which allow one to separate the S and P-wave components of the Dalitz plot. An 
asymmetry variable X^cp is then defined as 



X,CP - 1 2 ^ 1<5J 

with a defined as before. Legendre polynomial moments of orders greater than eight 
are neglected. A \ 2 f° r consistency with no CPV is calculated in the same way, this 
time accounting for correlations between the eight different angular moments I with 
a correlation coefficient obtained from simulation: 

x 2 = EEEW 3 (4) 

j=0j=0 

Here k is the number of intervals of invariant mass and the number of degrees of 
freedom is 8k. BaBar find the one-sided confidence level for no CPV to be 28.2% 
for the TT+vr", 28.4% for the vr+vr , 63.1% for the K+K, and 23.8% for the K+ir° 
sub-systems. It should be noted that this procedure assumes there is no interference 
from crossing channels and that only 5* and P waves are present. 



3 D + -> K~K + 7T + at LHCb 

A similar analysis strategy was pursued with the first 35 pb -1 of data taken at the 
LHCb experiment in 2010 [7]. The hadronic environment and the higher statistics 
present some additional challenges. For example, all multi-body analyses are subject 
to some extent to the variation of interaction asymmetries of particles with detector 
material across the Dalitz plot, because these asymmetries can depend on momentum. 
At LHCb the high boost tends to reduce these problems. This was demonstrated 
in the convenient Cabibbo-favoured control channel, — > K~K + 7T + , which has 
similar kinematics and signal yield to the D + decay. The reconstructed K~K + tt + 
mass distribution and D + — > K~K + tt + Dalitz plot is shown in Figure [2] There are 
around 370,000 D + — > K~K + it + signal events with 90% purity in the 2010 dataset. 

Monte Carlo simulations based on an amplitude model from CLEO-c [8] are used 
to develop binning schemes for the Dalitz plot with greater sensitivity to plausible 
types of CP violation than a simple regular grid. For example, it is found that the 
strong phase can change across resonances. If one bin is used for the whole resonance, 
the CPV would change sign within it and the Sep value for that bin would probably be 
consistent with 0. Therefore, without sacrificing model independence, the resonances 
are divided up into bins an appropriate way to account for the changing strong phase. 
With 25 bins in the Dalitz plot, the Monte Carlo studies suggest, for example, that a 
CP-violating phase in the 0(1020) of 5° would be observed at the 3cx C.L. with ~ 90% 
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Figure 2: K~K + n + Mass distribution (left) and D + -> K-K + n + Dalitz plot (right). 

probability. Other possible CPV signals are explored and, for example, the method 
is similarly sensitive to an ~ 11% CPV in the magnitude of the k(800). 

The charge asymmetries in the Df — > K~K + tt + control channel, and also in the 
Cabibbo-favoured D + — > K~h + 'k + decay, are investigated to eliminate the possibility 
of observing detector asymmetries as signals of CPV. No significant charge asymme- 
tries are observed, indicating that no Dalitz-plot dependent fake asymmetries are yet 
detectable. 

The signal mode is then binned according to four binning schemes, two of which 
account for the resonant structure of the decay and two of which, following the analysis 
at BaBar [5], do not. The Sep values in each bin are calculated and the resulting x 2 
for consistency with no CPV is converted into a p- value. The Sep values are shown 
in Figure [3] and the p- values are tabulated in Table [TJ No evidence of CPV is seen. 




Figure 3: Distributions of Sep variable for the LHCb D + — > K K + n + Dalitz plot, 
in phase space (top) and as pulls (bottom) for four choices of Dalitz plot binning. 
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Binning 


xViidf 


p-value (%) 


(a) Adaptive I 


32.0/24 


12.7 


(b) Adaptive II 


123.4/105 


10.6 


(c) Uniform I 


191.3/198 


82.1 


(d) Uniform II 


519.5/529 


60.5 



Table 1: Table showing x 2 /ndf and p- values for consistency with no CPV for the 
binnings shown in Figure [3] 



4 Four body decays 

The five-dimensional phase space of four body decays presents new challenges and 
new opportunities. The binned search for CP violation can be adapted to four-body 
decays, as done recently at LHCb in around 180,000 D° — > decays [9]. In 

four-body decays, the method of T-odd moments also becomes available. This has 
been applied in several analyses, for example by the BaBar Collaboration in their 
search for CPV in D + — )■ KgK + ir + 7r~ PE|. 

The proceedure is to define a triple product of the momenta of three of the par- 
ticles. 

C T = Pk+-(Ptt+ x (5) 
The fourth momentum is completely constrained. Then one defines an asymmetry 

At = r(c r >Q)-r(c T <o) 



r(c T > o) + v{c T < o) 



for D + decays and the analogous A^ for D~ decays. The quantity A = \{At — Aj) 
is then a measure of T violation. At and A^ can be individually nonzero if there are 
final state interactions [llj, but their difference cannot, unless CP is violated. In a 
sample of 20,000 decays, BaBar obtain A = (-1.2 ± 1.0 ± 0.46)%. 



5 Unbinned techniques 

There is a promising new technique available for unbinned model-independent searches 
which has not yet been applied to experimental data [12"] . It is possible to compare 
the densities of oppositely charged points in phase space using a two-sample test. A 
test statistic T is used to correlate the difference between D + and D~ p.d.f.s, f{x) 
and f(x'), 

T=\J J (f(x)-J(x))(f(x')-J(x'))^(\x-x'\)dxd^ (7) 
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To calculate T from data, one then translates this to: 

\ n 1 ™ 1 n '™ 

t = ^—^ j:m - 3D + ^— T ev(i^ - 3D - -= e ^ -^i) (8) 

where n (n) is the number of events (c.c. events). A sensible choice for weighting 
function ip(\xi — Xj\) is e~^ Xi ~ Xj ^ l 2a where the choice of a depends on the scale of 
the asymmetries expected in the Dalitz plot. 

The permutation test is then used to determine a p-value for consistency with 
no CPV. Many datasets where the "+" and "— " labels are assigned randomly to D 
decays are produced, such that there are n "+" decays and n "— " decays in each set. 
Then the p- value is the fraction of elements in the set for which T > T b s . 

Preliminary Monte Carlo studies indicate that this method should be more sensi- 
tive than binned searches, although it could also be very computationally expensive 
if applied to large datasets in its most basic form. 



6 Outlook 

Model-independent techniques for CPV searches are likely to become more important 
for the very large Belle II and LHCb data samples. Whilst they are not substitutes 
for full amplitude analyses, they provide the most straightforward means to establish 
CPV in a multi-body decay. It is hoped that model-independent techniques can shed 
light on some of the most pressing issues in charm physics in the near future. 
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